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Abstract. Multi— connected Universe models with space idenfication scales smaller 

^ ' than the size of the observable universe produce topological images in the catalogs of 

00 ' cosmic sources. In this review, we present the recent developments for the search of 

i the topology of the universe focusing on three dimensional methods. We present the 
crystallographic method, we give a new lower bound on the size of locally Euclidean 

1 multi— connected universe model of 3000 h' 1 Mpc based on this method and a quasar 

f*^ ■ catalog, we discuss its successes and failures, and the attemps to generalise it. We 

f™*^ ' finally introduce a new statistical method based on a collecting correlated pair (CCP) 

' technique. 
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1.1. Historical elements 



In the standard cosmological framework, the universe is decribed by a Friedmann- 
Lemaitre solution, the spatial sections of which being usually assumed to be simply- 
connected. Einstein equations being local, they allow us to determine the local 
geometry but they give no complete information about the global structure of the 
universe, i.e. about its topology, even if the geometry constrains to some extent 
the topology. This was pointed out by Friedmann himself just after he proposed 
its cosmological solution The indeterminacy about the global topology of 

the cosmological solutions was raised out as soon as Einstein proposed the first 
cosmological solution of his equations ; the Einstein static universe assumed spatial 
sections with the topology of the hypersphere S 3 , although de Sitter stressed that 
the same geometry could admit as well the projective space P 3 = S 3 /Z 2 as spatial 
sections ||. These two solutions were locally identical (same metric) but differed by 
their topology, i.e. by the choice of boundary conditions. 

Many arguments were then advanced in favour of a simply-connected universe, 
such as the simplicity or economy principle stating that one should introduce as 
few parameters as possible in physical modelling. Indeed, such an argument is very 
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unclear. There has been some tendency to favor spaces with finite volume. The 
eternal and infinite space of Newtonian physics was for instance leading to logical 
difficulties such as the Olbers paradox Q , and the finite cosmological solution proposed 
by Einstein in 1917 was received as a smart way of solving such paradoxes. The 
Mach's principle ^ |(| based on the idea that the local inertia is determined by the 
distribution of masses in the whole universe, also tends to favor universes with a finite 
volume. If space is simply-connected then it is finite if and only if it is locally elliptic, 
whereas if multi-connected it can also be locally hyperbolic or Euclidean. It has also 
been argued that infinite spaces were unacsthctic since "all phenomenon with a non 
vanishing probability must happen somewhere else" @, |) . This argument was used by 
Ellis to conclude that if space was infinite then one can avoid the former conclusion by 
dropping the assumption of spatial homogeneity and by arguing that we were living 
in a particular place of the universe [^). To finish, many arguments coming from 
some ideas in quantum cosmology are in favor of a finite volume space [H], [H], [l2| . 
Indeed, none of these arguments can help us to determine the shape and size of our 
universe. So, we are lead to the question, can we detect or constrain observationnally 
the topology of our universe ? 

Astonishingly enough, our century has on one hand seen the birth of the 
geometrical description of the universe and of a non static curved spacetime. On the 
other, mathematicians have developped the classification of 3D-manifolds. The first 
motivation came from crystallography, which lead to the classification of Euclidean 
3D-manifolds H [lj] and was achieved in 1934 |jl5). The classification of locally 
elliptic 3-manifold was set by F. Klein Jl6| and W. Killing [jl7| and solved by J. A. 
Wolf in 1960 |l8|. The classification of locally hyperbolic manifold was started in 
the 70's by W. Thurston (l|, |o) and is not yet achieved @ [22| f|§. Thus both 



elements necessary to answer the question of cosmic topology were developped at the 
same time without interacting so much. For more details concerning this interesting 
modern quest, one can see p4| JM p6L p7f. 



1.2. Mathematical elements 

In relativistic cosmology, our universe is described by a globally hyperbolic 4-manifold 
M. = £ x R [^8[ p9| , where the spatial sections S are homogeneous and isotropic 
Riemannian 3-manifolds. From a topological point of view, it is convenient to decribe 
such a manifold by its fundamental polyhedron (hereafter FP), which is convex. Its 
faces are associated by pairs through the elements of a holonomy group T which is 
acting freely and discontinuously on £ (see J3(], [5l], [3^] for mathematical definitions 
and J33|, |34|, |35| for an introduction to topology in the cosmological context). The 
holonomy group is isomorphic to the fundamental group 7Ti(£). Using the property 
(see e.g. f§): 

7Tl(M) = 7Tl(E X R) ~ m{R) X 7Ti(£) ~ 7Ti(£), (1) 

we deduce that the study of the topology of the universe reduces to the study of the 
topology of its spatial sections. 

What are the allowed homogeneous 3-manifolds usable in cosmology ? According 
to the sign of the spatial curvature K, the universal covering space of the spatial 
sections (which will be referred to as £) can be described by the Euclidean space E 3 , 
the hypersphere S 3 or the 3-hyperboloid H 3 if K vanishes, is positive or negative 



respectively. Thus, the homogeneous and isotropic 3-manifolds will be of the form 

E 3 /T, S 3 /T and H 3 /T. (2) 

Let us just summarize some of the properties of each family and also introduce G, 
the full isometry group of £ keeping the metric invariant. Indeed T will be a discrete 
sub-group of G. 

(i) locally Euclidean 3-manifolds: The covering space in the Euclidean space £ = E 3 
and the isometry group is G = R 3 x 5*0(3). The metric can be written under the 
form 

ds 2 =a 2 ( V ){-dr 1 2 + dx 2 +X 2 dn}, (3) 

where rj is the conformal time, a(rj) the scale factor, X the radial coordinate and 
dfl the unit solid angle. The generator of T are the identity, the translations, the 
reflexions and the hclicoidal motions. 

These transformations generate 18 different spaces jl8| among which 17 are multi- 
connected and correspond to the 17 crystallographic groups jl3|]. 10 spaces are 
compact and among them 6 are orientable. The description of these spaces can 
be found in e.g. Q. 

(ii) locally elliptic 3-manifolds : The covering space in the elliptic space = S 3 and 
the isometry group is G — SO (A). The metric can be written as 

ds 2 = a 2 (77) {-drf + d X 2 + sin 2 X dQ} ■ (4) 

With the curvature radius as unit length, the volume of the hypersphere is 

Vol(5 3 ) = f 4irsm 2 X d X = 2tt 2 . (5) 
Jo 

Wolf Jll| gives an exhaustive description of the allowed discrete groups T. They 
are the cyclic groups of order p {p > 2), the dihedral groups of order 2m, and the 
symmetry groups of the tetrahedron, octaedron and icosahedron. 

If we denote |T| the order of the holonomy group, it is straigthforward to show 
that 

Voi(s 3 /r) = 2^ 2 /|r|, ( 6) 

which tells us that the volume of an elliptic manifold is a topological invariant 
(degenerate since two different groups with same order will have same volume). 

(iii) locally hyperbolic 3-manifolds: The covering space in the hyperbolic space S = 
H 3 and the isometry group is G = PSL(2,C) = SL(2,C)/Z 2 . The metric can 
be written as 

ds 2 = a 2 (ry) {-drj 2 + d X 2 + sinh 2 X dtt} ■ (7) 

The classification of compact hyperbolic manifolds is not achieved yet and we 
will just describe two of them and give some useful properties. This classification 
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relies on the rigidity theorem |40|, [iT] stating that the geometry is fixed by the 
topology, a consequence of which being that the volume and other characteristic 
lengths are topological invariants, so that 



TiipO ~ m{Y) Vol(X) = Vol(Y). (8) 
It has also been shown that there is a minimal allowed volume [{42| 

Vol(E) > Vol min = 0.166. (9) 



We also define the outside radius r+, the radius of the smallest geodesic ball 
that contains the FP, the inside radius r_ , the radius of the biggest geodesic ball 
contained in the FP, and the injectivity radius rj n j, half the length of the shortest 
closed geodesic. r + , r_ and r- m j are altogether topological invariants. 

The classification and the description of the known compact hyperbolic manifolds 
can be obtained by using the software SnapPea |3§] which provides the FP, the 
generators of the holomy group and the characteristic lengths. Examples of such 
manifolds can be found in e.g. [^| . Here we mention for a later use a manifold 
of special interest, the Weeks space [p7| , which is the smallest known compact 
hypermabolic manifold. Its FP has 18 faces and its geometrical characteristics 
are 

Vol = 0.94272, r + = 0.7525, r_ = 0.5192, r inj = 0.2923 (10) 



1.3. Physical and observational introduction 

At the time being, two classes of methods to detect and/or constrain the topology of 
the spatial sections of our universe have been proposed and investigated. They use 
two classes of cosmological observations, namely the 2D cosmic microwave background 
data and the 3D catalogs of discrete sources. 

(i) 2D-mcthods : The cosmic microwave background is composed of all the photons 
emitted during the decoupling period between matter and radiation at a redshift 
of z ~ 1100. It is observed as a black body with a temperature of 2.7 K with 
small fluctuations of order 10~ 5 . The temperature anisotropies reflect the small 
inhomogeneities that ultimately lead to the observed structures in the universe. 

When studying these anisotropies in a multi-connected universe, one has to take 
into account the fact that we are in a compact manifold and thus that the 
wavelengths must be discretized. Stevens et al. |i"7| showed that in a cubic 
hypertorus one should see a cut-off in the two-point correlation function if 
L/Rh < 0.8. This analysis was then generalised to non cubic hypertorus 
and to the six flat manifolds ^0|. However, in compact hyperbolic manifolds 
there are super-curvature modes []l9| , which implies that no such cut-off exists 

nil- 

It has then been realised that the patterns on such an extended surface must 
be correlated if the universe is multi-connected. Levin et al. |5l[] studied these 
correlations in the flat manifolds and Cornish et al. [[l9| developped a topology 
independent method to look for a topological signature in the cosmic microwave 
background map as expected from the MAP and Planck satellites. 
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All these studies have been performed assuming that the small initial 
inhomogeneities where generated during an inflationary phase. There is however 
another process that could generate such inhomogeneities, namely topological 
defects produced during a symmetry breaking phase transition in the early 
universe. Uzan and Peter ]54| showed that the topology implies a constraint 
on the topological defects network from which it was deduced |55| that there will 
be a cut-off in the angular power spectrum of the temperature anisotropies, the 
cut-off value depending only on the characteristic size of the universe and on 
the cosmological parameters. Note that such considerations also provide another 
method for constraining the topology, namely the observation of an extended 
defect at a redshift z will give a lower bound on the size of the universe. 

In conclusion, at the moment 2D-methods constrain flat manifolds to L/Rh < 
0.8, and there is no convincing constraint on compact hyperbolic manifolds |56|| . 
Some powerful methods can be hoped to be used when high resolution cosmic 
microwave background maps are achieved. 

(ii) 3D-methods : In the framework of pre-relativistic cosmology, it was already 
pointed out by Schwarzschild that if space is multi-connected, one could see 
multiple images of astronomical objects p3j , whereas Friedmann jO did the same 
in the framework of general relativity. All 3D-methods for testing the space 
topology are based on this fact. We shall develop their description in the following 
sections. 

In these Proceedings, we describe the developments of the 3D-mcthods, from 
the Schwarzschild initial idea of the existence of multiple images of the same object, 
until its statistical implementation under the form of cosmic cristallography. We then 
carefully study the method and explain why, contrary to what was previously thought, 
it won't work in locally hyperbolic universes. We then present attempts to generalise it 
and finish by describing a new and promising "CCP" method. This will be examplified 
by numerical simulations with depiction of color pictures in order to make the subject 
as clear as possible. 



2. FROM THE ORIGINAL IDEA TO THE CRYSTALLOGRAPHIC 
METHOD 

2.1. Initial idea 

As stressed above, if space is multi-connected, one should see multiple "topological" 
images of the same object. This is illustrated on figure [l] where we have used a 2-D 
example. In the universal covering space, the spacetime trajectories of the different 
images (A, B, . . .) of a given object will intersect the observer past light cone at 
different times (a, b,...). The same object will then been seen at different stages 
of its evolution. Indeed in a 4D universe, these images will also be seen in different 
directions of the sky. 

The question is indeed to find methods to implement such a property in order to 
detect or to put lower bounds on the size of the universe. 
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Figure 1. a 1 + 1 analogy of an multy-connected universe S 1 X R. In the universal 
covering space, we have plotted the past light cone of the observer (blue), the 
trajectories (green) of different topological images (A-F) of a given object, assuming 
it is comoving, and the spacetime point (a-f) where they are observed. O will see 
them at different look-back times. Indeed in 3 + 1 dimension, he will also see them 
in different directions in the sky. 



2.2. Multiple images of individual objects 

Some specific objects have been tried to be recognised individually. An extensive 
description of these attempts can be found in J33| and we just stick to a brief 
description. 

• Milky Way: Can we recognise our own Galaxy ? The maximum distance up to 
which we would be able to do that has been estimated as 7.5 Mpc by Sokolov 
and Shvartsmann [j58| which gave the observational bound r + > 15 /i -1 Mpc. 
Then, Fagundes and Wichoski |HJ proposed that our galaxy was a quasar long 
ago but indeed, for physical reasons due to realistic quasar models, our Galaxy 
cannot have been a quasar in its past. Fagundes [jfl) also studied the occurrence 
of images of the Milky Way in the Best hyperbolic space model J74| . Presently, 
no source has been identified as an image of our own Galaxy. 

• Galaxy clusters : Gott used the Coma cluster [6^] has a candidate for the search 
of ghosts images and obtained the bound r + > 60 /i -1 Mpc. He also performed 
simulations in T 3 to provide a pattern of clustering and a correlation function in 
agreement with observations of nearby galaxies. Roukema and Edge |4(| proposed 
3 candidates as ghost images for the Coma cluster. Further investigations are 
needed. 

Sokoloff and Shvartsmann |58| considered the Abell catalog of clusters and gave 
the bound r + > 600 hr 1 Mpc in Euclidean space, whereas Demianski and Lapucha 
J60| searched unsuccessfully for opposite pairs in a catalog of 1889 clusters. 
Lehoucq et al. developped a statistical method that is presented below to obtain 
r + > 650 /i -1 Mpc in Euclidean space and Fetisova et al. reported a spike about 
125 hr 1 Mpc in the correlation function of rich clusters, but such an evidence does 
not necessarily supports multi-connectedness (see e.g. [fj3|). 

• Quasars: Quasars occupy a large volume of the universe [z ~ 3). However, they 
are probably short-lived compared to the expected time necessary for a photon 
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to go around the universe. If T q ~ 10 9 years, then one can hope to test topology 
on scale ~ r + ~ cr^. So, they are of interest even if r q ~ 10 8 years. Paal [ pEfl 
remarked that although individual quasars may have short lifetime, they may be 
part of larger associations which survive longer. 

Roukema p5[ tried to identify quintuplets and quadruplets of quasars in a catalog 
of 4554 quasars. He found 27 identifications whereas he found 26±1 identifications 
in simulations of a simply-connected universe. Fagundes |7lj tried to discuss the 
redshift contreversy in the Best |74j model. 

Until now, the best bound obtained by 3D-methods was r+ > 650/i _1 Mpc and 
hold only in Euclidean universe. There is still some observational room for the 
topologies of an hyperbolic universe, even on scales significantly shorter than the 
horizon radius. 

All the methods trying to recognize a particular object suffer from the same major 
problems: 

(i) distance determination: catalogs give the angular position and the redshift. One 
has then to go from redshift space to real space, which depends on the knowledge 
of the cosmological constant and of the density parameter. This explained why 
most of the former studies were performed assuming tto = 1, J7a = 0. This point 
will be develop later. 

(ii) peculiar velocities: objects are not comoving, so that the position of their 
topological images are shifted. If their velocity is of order v p ~ 500km.s , then, 
assuming that the time for a photon to go around the universe can be estimated 
by t ~ r+/c, one obtains the maximal precision on the determination of position 

- (vp/c) x r+. 

(iii) morphology: when one tries to associate two images to the same physical object, 
one has to be aware that, in general, the object will not be seen in the same angle 
and at the same stage of its evolution. It will then be very difficult to recognise 
it and the identification will strongly depend on galactic evolution models. 

Conversely, the discovery of topology by any means could provide informations 
such as peculiar velocities with a better accuracy ]7q] . We will also be able to see 
the same object at different stages of its evolution and on different angles, which will 
constrain galactic evolution models. 



2.3. The crystallographic method 

Since there is little chance to recognise different images of a given object, one can 
try to detect these images statistically, since for instance the number of topological 
images of a single object in a toroidal universe has been estimated in table [l] of [Q. 

The crystallographic method is based on a property of multi-connected 
universes according to which each topological image of a given object is linked to 
each other one by the holonomies of space. Indeed, we do not know these holomics 
as far as we have not determined the topology, but we know that they are isometries. 
For instance in locally Euclidean universes, to each holonomy is associated a distance 
A, equal to the length of the translation by which the fundamental domain is moved to 
produce the tessellation in the covering space. Assuming the FP contains N objects 
(e.g. galaxy clusters), if we calculate the mutual 3D-distances between every pair of 
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Number of topological images of a single object 



L 



N(z < 1000) 



N(z < 4) 



N(z < 1) 



500 



7000 



1200 



180 
7 



1500 



279 



45 



2500 



GO 



10 



1.5 



Table 1. Number of topological images of a single object in a toroida 
characteristic size L (in h —1 Mpc). Reprinted from Lehoucq et al. 5- 



universe with 



topological images (inside the particle horizon), the distances A will occur N times 
for each copy of the fundamental domain, and all other distances will be spread in a 
smooth way between zero and two times the horizon distance. In a histogram plotting 
the number of pairs versus their 3D separations, the distances A will thus produce 
spikes. Simulations indeed showed that the pairs between two topological images of 
the same object drastically emerge from ordinary pairs |34| in the histogram. The 
applicability of the method in Euclidean spaces has also been discussed by Fagundes 
and Gaussmann f33f when the size of the physical space is comparable to the horizon 
size. 

Two kind of catalogs of astronomical objects can be thought of to apply this 
method: the galaxy cluster catalogs, which typically have a redshift depth z = 1, and 
the quasars catalogs, which typically extend to z = 3. Concerning quasars, even if 
their lifetime is probably too short to be good candidates for producing topological 
images, they are usually part of systems that have a much larger lifetime [[35]. Let us 
stress that a recent survey |36| in the Hubble Deep Field south NICMOS field found 
17 galaxies with a redshift between 5 and 10 and 5 galaxies with a redshift above 10, 
among a total of 323 galaxies. This can let us hope to apply this method to deeper 
catalogs in the future. The angular resolution needed is given by the fact that the 
objects have a peculiar velocity and that they will not be seen at exactly the same 
position Q. Note that the crystallographic method, contrary to the "direct" method 
which would try to recognize topological images of individual objects, is not plagued 
by the fact that topological images of the same object are seen at different stages of 
its evolution. 

We present on figure || an example from p3 for the flat torus with a cubic FP. 

Now, we apply this method to a catalog |7q] containing more than 11,000 quasars 
up to a redshift of z ~ 3.25 (see figure || where we have depicted its projection on 
the celestial sphere and the redshift distribution of its objects). The pair separation 
histogram is plotted on figure |] and it exhibits no topological signature. Since 
respectively 60 %, 80 % and 95 % of the sources are below the redshifts 2,2.3 and 3, 
we only test scales smaller that z ~ 3. 

In a locally Euclidean universe, the distance-redshift relation is given by [77| 



We simulated a catalog with the same number of objects and the same redshift 
distribution to reproduce the real catalog depicted in figure 0. Varying the size of 
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Figure 2. Application of the crystallographic method in a cubic torus universe, 
from [Q. Here, A corresponds to the separation distance and V to the volume of 
the manifold. 



cubic hypertorus, a spike appears when Lq ~ 3000 h 1 Mpc for the simulated catalog 
and thus we obtain the new bound on the size of a locally Euclidean manifold as 

L a > 3000 h~ l Mpc, (12) 

which is the best constraint presently available with 3D-methods for locally Euclidean 
manifolds. Note that it corresponds to 



Rh 



> 0.5, 



(13) 



which is comparable to the bound obtained with 2D-methods. Note that this result 
has not been published elsewhere. 
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Figure 3. The quasar catalog: (left) the distribution of objects on the celestial 
sphere in Hammer-Aitoff equal area projection and (right) the redshift distribution. 



3. WHY THE CRYSTALLOGRAPHIC METHOD DOES NOT WORK 
IN HYPERBOLIC SPACES 

When one looks carefully at the crystallographic method, one must wonder about the 
origin of the spikes. As it has recently been explained pffl, two kinds of pairs can 
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3D-dis1ance 

Figure 4. The pair separation histogram for the quasar catalog (see figure^) exhibits 
no topological signature when assuming a locally Euclidean geometry. 

create a spike, namely (we keep the notation and vocabulary introduced in p7|) 

(i) Type I pairs of the form {g{x), g(y)}, since dist[g(ir), g(y)] — dist[x,y] for all 
points and all elements g of T. 

(ii) Type II pairs of the form {x,g(x)} if dist[a;, g{x)\ = distky, g{y)\ for at least some 
points and elements g of T. 

Both families of pairs are depicted on figure [| for the example of the 2-torus. 
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Figure 5. Difference between type I and type II pairs in the example of the two 
dimensional torus. The translations gi, 32 Sj -1 and g^ 1 are the generators of the 
holonomy group. Type I pairs are the ones between the ghosts of two distinct objects 
(hexagon or square), and type II pairs are the ones between two topological images 
of the same object. 



3.1. Clifford translations 

A Clifford translation is an isometry g such that the displacement function dist(ir, g(x)) 
is constant Jigf ] . This is precisely what is required to get type II pairs in the histogram. 

In compact hyperbolic manifolds, dist[x, g(x)] always depends on the position x of 
the source (see [67 for details). 

In elliptic spaces, distances are position-independent whenever the holonomy is a 
Clifford translation. All finite groups of Clifford translations of spheres are the cyclic 



group, the binary dihedral, tetrahedral, octahedral and icosahedral groups |69 , from 
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which we deduce 67 that the covering transformations which move a source to its 
nearest neighbours are Clifford translations (although the transformations to more 
distant neighbours might not be), and type II pairs can be produced. 

In locally Euclidean universes, type I and type II pairs are both present. The 
reason is that the 3-torus has the very special property that the separation distance 
of gg-pairs (i.e. any pair of images comprising an original and one of its ghosts, or two 
ghosts of the same object) is independent of the location of the source. In other Eu- 
clidean spaces the spectrum of gg-pair distances varies with the location of the source. 
However all closed Euclidean 3-manifolds have the 3-torus as a covering space, so for 
each such manifold there will be some distances which are independent of the loca- 
tion of the source. As a consequence, the topological signal expected in the histogram 
from type I and type II pairs clearly stands out, as was shown in the simulations of f64j . 

In conclusion, type II pairs exist only in elliptic and Euclidean universes. The 
question we still have to answer to judge the efficiency of the crystallographic method 
in the general case is: what about the amplitude of the spikes related to type I pairs 
only ? 



3.2. Cosmological parameters 

If the injectivity radius of the space manifold is smaller than the calatog's length scale, 
type I pairs will always be present, whatever the curvature. Their number roughly 
equals the number of copies of the fundamental domain within the catalog's limits. 
Following |67f , this number can easily been estimated by computing the ratio between 
the volume of the geodesic sphere of radius x( z ) an d the volume of the manifold. The 
calculation (in the case of hyperbolic universes) leads to 

N(Q O -^7\ ^(^h2 X (Z)-2 X (Z)) 

^v(i2 m0 , s2a , z < zj — vbiTsj ' 

where x(^) is the radial distance corresponding to the redshift Z , given 



X(z) = / — = / ; = , (15) 



a aa 



TTT xJCIaoX 2 + (1 - fi m - ^Ao) 



where f2 m o and Q\ are respectively the density parameter and the cosmological 
constant. We have plotted on figure || an estimation of the number of copies of an 
object taken respectively in a catalog of galaxy clusters and in a catalog of quasars, 
as a function of the cosmological parameters. 

The small number of type I pairs in hyperbolic manifolds is due to the property 
that the volume of the manifold is fixed once the topology is determined (the rigidity 
theorem), contrarily to Euclidean spaces where the characteristic sizes and the volume 
of the FP can be chosen at will (since K = 0, the geometry does notimpose any 
characteristic size). 

In conclusion, type I pairs will indeed contribute to spikes in the pair separation 
histogram but the amplitude of these spikes will be of the same order than the 
statistical noise. 
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Figure 6. The number of copies of an object belonging to a catalog of clusters (left) 
and a catalog of quasars (right) as a function of the cosmological parameters. 



3.3. Numerical simulations 

The former conclusion can be tested numerically. For that purpose, we concentrate on 
locally hyperbolic three-dimensional manifolds. It can be first embedded in a four- 
dimensional Minkowski space by introducing the set of coordinates (x M ) M= o..3 related 
to the intrinsic coordinates (x, 9, p) through (see e.g. |6£], |70[ |) 

xq = coshx 

x\ = sinh x sin 9 sin tp 

X2 = sinh x sin 9 cos tp 

x 3 = sinh x cos 9, (16) 
so that the three-dimensional hyperboloid H 3 has the equation 
; +xl +xl = -1. 



r 2 

Xq 



(17) 

With these notations, the comoving spatial distance between two points of comoving 
coordinates x and y can be computed directly in the Minkowski space by [fri] 



d[x,y] = argcosh [x^y^] 



(18) 



where x M = r/^x" , r/^ being the Minkowskian metric coefficients. 

To generate an idealised catalog, C, we first distribute homogeneously N objects 
in the FP, then we by unfold the catalog by applying the generators of the holonomy 
group. Now, to generate a catalog with the required depth z max in redshift, the set of 
all points obtained as above must be truncated in order to keep only the points x such 
that dist[0, x] < xi^mo, ^ao! z max), given by equation (15). This amounts to select 
objects located within the geodesic ball of radius x(^mu> ^ao! z max ) centered onto an 
observer right at the centre. 

The unfolding of Weeks manifold is shown in figure 0, where we have depicted the 
central cell and its 18 nearest neighbours in Klein coordinates. In figure ^, we have 
plotted a simulated catalog showing all points, including those outside the horizon 
distance. 

We then compute all the 3D separations between all the pairs of C(z) and we plot 
the histogram of the number of pairs with a given separation. 
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Figure 7. Unfolding the Weeks manifold: we have represented the FP and its 18 
nearest copies in Klein coordinates. 
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Figure 8. Building a catalog C for the Weeks manifold. The green dots are objects 
located within the horizon. The red ones lie outside and will be rejected from C. 

Indeed the former procedure applies when the observer stands at the center of 
the polyhedron (x = 0). In |67|] , we explain how to take into account an off-centered 
observer. 

In Lehoucq et al. [ |67[ , we have generated some pair histograms in a universe whose 
spatial sections have the topology of the Weeks manifold, using N = 1000 objects in 
the FP. In figure || we give an example wher f^A = 0, f2 mu = 0.2 and where the observer 
stands at the center of the polyhedron (x = 0). The dependence on the cosmological 
parameters, on the position of the observer and on the catalog depth where studied 
in details in [ |67[ . 
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Figure 9. Pair separation histogram for a universe with the topology of the Weeks 
manifold, Co = 0.2, Ha = 0, for a simulated catalog of cosmic objects with depth 
Zmax = 1- No spike stands out. 



As expected, none of the plots exhibit spikes. This is due to the combined effects 
of the geometry and of the cosmological parameters values. 

To summarize the failure of cosmic crystallography for hyperbolic manifolds: 

(i) hyperbolic manifolds are such that dist[a;, g(x)] depends on x, so that there is no 
amplification for the type II pairs {x,g(x)}, whereas dist[;c, g{x)] = dist[y, g(y)} 
in the Euclidean case. This suppresses the corresponding spikes. 

(ii) The spikes associated to the isometries (i.e. such that V.g € T, dist[g (x), g(y)] = 
dist[x,y]) must remain. But, given the cosmological parameters, we have shown 
that the number of topological images is too low to create such spikes associated 
to type I pairs. 

4. GENERALISING THE CRYSTALLOGRAPHIC METHOD 

The next question is: can we generalise the crystallographic method in order to make 
it work in any case ? 

Two solutions were proposed almost immediately to improve the crystallographic 
method and to extract from 3D catalogs the signature of type I pairs: 

(i) Fagundes and Gaussman [f72[ substracted the pair separation histogram for a 
simulated catalog (with the same number of objects and the same cosmological 
parameters) in a simply-connected universe to the observed pair separation 
histogram. The result is "a plot with much oscillation on small scale, modulated 
by a long wavelength quasi-sinusoidal pattern" . The statistical relevance of this 
signature still has to be investigated. To illustrate this point we show on figure [lO] 
the result of a simulation where we have substracted to realisations on a simply- 
connected Friedmann-Lemaitre universe. This has to be compared with figure 1 
from [ [72[ and let us think that is method is not efficient for determining a relevant 
signature of the topology of the universe. 

(ii) Gomero et al. |6^] proposed to split the catalog in "smaller" catalogs and 
average the pair separation histograms built from each sub-catalog to reduce the 
statistical noise and extract a visible signal of the non-translational isometries. 
The feasability of this method has not yet been demonstrated. 
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Figure 10. The plot for the Weeks manifold [grey] compared to the difference 
between two realisations in simply-connected universes with same local geometry 
[green] . 

5. A NEW EFFICIENT STATISTICAL METHOD 

We now present a new method Q based on the property that whatever the topology, 
type I pairs will always exist. Indeed, as shown ahead, this does not lead to any 
observable spike in the pair separation histogram. We thus had to find a method 
to enhance the topological signal by "collecting" all the homologous pairs together. 
This can be achieved by building a collecting correlated pairs method (hereafter CCP). 
Such an approach (described in details in |73|] and below) is not able to determine the 
exact topology, but will provide a signature of the compactness of the spatial sections, 
which is indeed a first step toward the determination (or the rejection) of the cosmic 
topology. 

In this section, we recall the construction of the CCP-index and the main resuts 
we have obtained. 



5.1. Basic idea and general method 

As stressed in the previous section, type I pairs will always exist in multi-connected 
spaces as soon as one of the characteristic scales of the fundamental domain is smaller 
than the Hubble scale. Defining ff»li<j<2jc as the IK generators of T and referring to 
x as the position of the image in the universal covering space X, we have: 

(i) Vx,yeX,Vger, 

dist[g(x),g(y)]=dist[x,y\. (19) 
We will refer to these pairs as xy-pairs. 

(ii) Vx e X, Vgi,g 2 G I\ 

dist[gi(x),gi o g 2 (x)] = dist[x, g 2 (x)]. (20) 



We will refer to these pairs as x<7(a;)-pairs. 
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Both the xy-pairs and the a;g(x)-pairs are type I pairs. 

To collect all these pairs together and enhance the topological signal, we define 
the CCP-index of a catalog containing N objects as follows. 

(i) We compute all the 3D-distances dist[x, y] for all points within the catalog's limit. 

(ii) We order all these distances in a list di\ 1<i<p , where P = N(N — l)/2 is the 
number of pairs, such that di + \ > di. 

(iii) We create a new list of increments defined by 

Vie [1...P-1], \ = d l+1 -d, (21) 

(keeping all the equal distances, if any, in the list). 

(iv) We then define the CCP-index TZ as 

T^p^T, (22) 

where TV = card({i, A* = 0}), so that < K < 1. 

With such a procedure, all type I pairs will contribute to TV. For instance, if a given 
distance appears 4 times in the list di\ 1<i<Pl it will contribute to 3 counts in J\f. 
Compared to the old crystallographic method, all the correlated pairs are gathered 
into a single spike, instead of being smoothed out into the noise of the histogram pair 
separation. 

Indeed, in a more realistic situation, one has to take into account bins of finite 
width e and replace TV by 

W e = card({i, A, E [0,e[}) (23) 

in the computation of TZ. The effect of this "binning" is discussed below. We now 
focus on the "idealised" version of the procedure by studying the amplitude of the 
CCP-index in some multi-connected models. 

For that purpose, let us assume that the catalog is obtained from an initial set of 
A objects lying in the fundamental domain and that B copies of the domain are within 
the catalog's limits (B = if the whole observable universe up to the catalog's limit is 
included inside a fundamental domain). The total number of images is N = A(B + 1). 
Indeed B is usually not an integer but we assume it is, in order to estimate the 
amplitude of TZ and compare it with the result in a simply-connected model. In f73[| , 
it has been shown that 



A/'min — A 



(A-1)|+M(£,B) 



(24) 



where ^i(S, B) is a function characterising the manifold S. It is an increasing function 
of B which vanishes for B = 0. 

Indeed, if the holonomy group T contains Clifford translations allowing for type II 
pairs, or if there are "fake" pairs (i.e. such that dist[x,y] = distkt, «]), /V m i n computed 
from (|24|) will give a lower bound for the true Af. 

The normalised CCP-index (|2^) follows straightforwardly. TZ is a good index for 
extracting the topological signal since 
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(i) when B = (i.e. when the fundamental domain is greater than the catalog's 
spatial scale), 1Z = 0, 

(ii) when the number of sources in the fundamental domain becomes large, it behaves 
as 

R ^™^) m ^4 » co. p.) 

As shown in |f73f , one can compute analytically the CCP-index in various examples 
and show that it is a fair indicator of the existence of at least one compact spatial 
dimension. 



5.2. Simulations and statistical relevance 

Now that the CCP-index can be computed for any multi-connected space, we can 
compare it with its value in a simply-connected Friedmann-Lemaitre model containing 
the same number objects (i.e. A(B + 1)). 1Z is given as a function of V\, A (the number 
of sources in the fundamental domain) and B (the number of copies of the domain 
within the observable universe) by: 

K = m^-B)A-B\ 

(B + l) 2 A 2 -(B + 1)4- 2" 1 ' 

Indeed, when working with a real catalog, we do not know the radial distance of 
cosmic objects but only their redshift. As seen before, the determination of the radial 
distance requires the knowledge of the cosmological parameters Qq and A and can be 
obtained analytically only when A = 0. 

Thus, if the universe is multi-connected on sub-horizon scale, the plot of TZ in 
terms of fio and will exhibit a spike only when the cosmological parameters have 
the right value (as shown on plot 7 of |7^| ). If the cosmological parameters are not 



exactly known, the distance determination (15) will be wrong and the topological 



signature will be destroyed (see figure 6 in J73|). 
Two consequences follow: 

(i) One should span the parameters space (CIq, Qa) m order to detect the topological 
signal, plotting 7^.(f2o,^A)- 

(ii) If there is any topological signal, the position of the spike gives the values of the 
cosmological parameters on the scale of the catalog's limit (see figure |TT| ) . 

We proceed as follows. We first generate a catalog by choosing the number of 
objects in the fundamental domain (A — 30), the topology (Weeks manifold) and the 
cosmological parameters (e.g. £Iq = 0.2, £1a = 0.1), we then use a second code to 
apply the test and we draw 1Z in terms of the two cosmological parameters. As shown 
in jf3|, the method works pretty well in the sense that there is a strong spike which 
signals a non trivial topology and which determines the cosmological parameters. But 
we also see that a slight deviation in the evaluation of the cosmological parameters 
would make the spike to disappear. This effect is now discussed. 
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5.3. Real data 

When one wants to apply the CCP-method to real data, one has to face a number of 
problems. 

First, we cannot use a zero width bin, one of the reasons being that the sources 
are not comoving. In |f73| , we estimate the precision needed on the cosmological 
parameters when working with a bin width e such that 



sn 



(27) 



Indeed, using a catalog with a smaller depth z max = 3 will allow us to use a smaller 
resolution for the cosmological parameters. One has thus to find a compromise between 
depth and resolution as discussed in |7^] . On figure O, we give an example with a bin 
width of e — 1CP 6 which produces a background noise. 

Now, we apply our test to the quasars catalog 179] used ahead (see figure 0). 



No topological signature was found 73 . Does it mean that there is no topological 
effect on scales smaller than z max ~ 3 ? Not necessarily, since we applied the test 
with presicions e = 10~ 7 , 1CP 6 , 10~ 5 and were unable to span the full cosmological 
parameters' space with the required accuracy. The computational time is one of the 
main limitations of our technique (this point is discussed extensively in J73[) 

R-histogram (e=10~ 6 ) 



30 
25 




Figure 11. Computation of the CCP— index on a simultated catalog of depth z = 3 
in a universe model with the Weeks topology, Qo = 0.2, Qa = 0.1, using a bin 
resolution e = 10 — 6 . 

Now, another limitation comes from the peculiar velocities of the sources. In ]73[| , 
we discussed in details the effect of peculiar velocities both on icy-pairs and on xg(x)- 
pairs, we and showed that even if the amplitude of the topological signal is reduced, 
it is not destroyed. 

6. CONCLUSIONS 

We have reviewed the main allowed topological structures for relativistic universe 
models and we discussed the various methods aimed to detect this topology. Focusing 
on 3D-methods, we described the crystallographic method and we applied it to 
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a quasars catalog to obtain a new constraint for Euclidean manifolds, L > 
3000 h^ 1 Mpc. Then we discussed the failure of this method for detecting the topology 
of compact hyperbolic manifolds. We thus generalised the method by introducing a 
new technique based on the construction of a CCP-index. The main difference with 
the cosmic crystallography method is the collecting process of all correlated pairs of 
the catalog, which enhances the signal associated to the existence of a non trivial 
topology. We gave examples of computation of such an index in order to show its 
statistical relevance. We finally discussed the implementation difficulties of this new 
method when working with observational data. 
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